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1.  INTRODUCTION 


1.1  Project  Objectives 

The  objectives  of  this  two-year  study  are: 

(1)  Assemble  three  data  sets  to  be  used  to  test  and  evaluate  the  performance  of 
neural  networks  for  automated  processing  and  interpretation  of  seismic  data 
(Table  1). 


Table  1 .  Seismic  Data  Sets  for  the  DARPA  Neural  Network  Program 


Data  Set 

Description 

Data  Set  #7 

Data  from  approximately  300  events  to  develop  and  train 
neural  networks  to  perform  seismic  data  processing  and 
interpretation  tasks  such  as  automated  phase  association, 
onset  time  estimation,  typical  and  atypical  event  recogni¬ 
tion,  and  event  identification  [LaCoss  et  al.,  1990]. 

Data  Set  #2 

Data  from  approximately  30  events  to  test  the  response  of 
the  neural  networks  to  novelty  signals.  These  data  are 
recorded  at  the  same  stations  as  the  events  in  Data  Set 
#1,  but  are  from  different  source  types. 

Data  Set  #5 

Data  from  approximately  300  events  to  test  the  generality 
and  adaptability  of  the  neural  networks.  These  events  are 
recorded  by  stations  in  a  different  geologic  environment 
than  the  stations  used  for  Data  Set  #1. 

These  data  sets  arc  to  be  provided  to  a  group  at  MIT  Lincoln  Laboratory  who  is 
developing  and  testing  neural  networks  for  the  reismic  application  of  the  DARPA 
Neural  Network  Program. 

(2)  Evaluate  the  results  of  the  neural  network  program  in  the  context  of  monitoring 
nuclear  explosion  testing. 

1.2  Current  Status 

Much  of  our  effort  during  the  first  year  of  this  project  was  on  the  development  of 
Data  Set  #1.  This  data  set  consists  of  short-period  waveforms  and  parametric  data 
from  241  regional  events  recorded  by  the  NORESS  and  ARCESS  arrays  in  Norway. 
Data  Set  #1  also  includes  parametric  data  from  249  other  events  (e.g.,  arrival  times, 
amolitudes,  polarization  attributes,  etc)  recorded  at  NORESS  and  ARCESS  during  a 
continuous  10-day  period.  The  delivery  of  this  data  set  to  MIT  Lincoln  Laboratory 
was  completed  in  March,  1991  (the  total  data  volume  is  about  1.2  GBytes).  Since 
their  project  began  several  months  before  ours,  we  also  provided  them  with  3- 
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component  waveform  data  recorded  by  the  center  element  of  the  NORESS  array  from 
73  regional  events  that  we  assembled  under  a  separate  DARPA  contract  These  data 
were  delivered  to  MIT  Lincoln  Laboratory  in  March  1990. 

We  recently  began  to  assemble  Data  Set  #2  which  will  consist  of  about  30  earth¬ 
quakes  recorded  at  NORESS  and  ARCESS  (only  18  of  the  241  events  in  Data  Set  #1 
are  earthquakes).  We  expect  to  deliver  this  test  data  set  to  MIT  Lincoln  Laboratory  in 
September,  1991.  Data  Set  #3  will  probably  consist  of  events  recorded  by  regional 
arrays  in  Germany,  GERESS,  and  Finland,  F1NESA.  Data  from  these  areays  are 
currently  being  archived  at  the  Center  for  Seismic  Studies  (CSS).  We  will  assemble 
this  data  set  during  the  second  year  of  our  project. 

We  have  several  concurrent  efforts  directed  towards  the  evaluation  of  neural  net¬ 
work  techniques  in  the  context  of  monitoring  nuclear  explosion  testing.  First,  we  plan 
to  integrate  the  neural  network  developed  by  MIT  Lincoln  Laboratory  for  automated 
regional  phase  association  into  the  Intelligent  Monitoring  System  (IMS).  IMS  is  a 
DARPA-sponsored  computer  system  for  automated  processing  and  interpretation  of 
seismic  data  recorded  by  arrays  and  single  stations  [Bache  et  al.,  1990,  1991].  This 
system  has  been  in  operation  at  CSS  since  October,  1990  (its  predecessor,  the  Intelli¬ 
gent  Array  System,  has  been  in  operation  since  October,  1989).  The  neural  network 
developed  by  MIT  Lincoln  Laboratory  assigns  a  regional  phase  identification  Kg.,  Pn , 
Pg,  Sn,  Lg,  or  Rg)  to  detections  registered  at  array  stations.  We  will  integrate  this 
neural  network  into  IMS,  and  compare  its  performance  to  that  of  the  current  rule-based 
expert  system  using  data  recorded  at  NORESS  and  ARCESS.  MIT  Lincoln  Laboratory 
sent  us  their  software  module  in  August,  1991.  We  are  currently  testing  this  module, 
and  we  will  begin  system  integration  within  the  next  few  weeks. 

Another  important  problem  in  automated  seismic  data  interpretation  is  initial 
phase  identification  (P  or  S)  using  data  recorded  by  3-component  stations.  We 
developed  neural  networks  for  this  application,  and  trained  and  tested  them  on  data 
recorded  by  the  3-component  elements  of  the  NORESS  and  ARCESS  arrays,  and  on 
data  recorded  by  the  3-component  IRIS  stations  in  the  Soviet  Union.  This  effort  is 
described  in  detail  in  Volume  II  of  this  report  [Patnaik  and  Sereno,  1991b].  We  are 
integrating  this  neural  network  into  IMS,  and  we  will  compare  its  performance  to  that 
of  the  rule-based  system  using  data  recorded  by  the  IRIS  stations. 

1.3  Outline  of  the  Report 

This  annual  report  is  divided  into  two  separate  volumes.  Volume  I  (this  docu¬ 
ment)  is  a  description  of  Data  Set  #1  that  was  provided  to  MIT  Lincoln  Laboratory  for 
their  neural  network  study.  Volume  II  presents  the  results  of  our  own  neural  network 
application  to  the  problem  of  initial  phase  identification  using  polarization  attributes 
derived  from  3-component  data  f Patnaik  and  Sereno,  1991  b). 

The  main  technical  section  of  Volume  I  is  Section  2.  Section  2.1  describes  the 
NORESS  and  ARCESS  arrays  and  instrumentation.  Section  2.2  gives  a  description  of 
the  regional  events  in  Data  Set  #1  (e  g.,  location,  magnitudes,  distances,  and 
identification).  The  exchange  format  for  waveform  and  parametric  data  is  described  in 
Section  2.3  and  Appendix  A.  Section  3  summarizes  Data  Set  #1. 
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2.  DATASET#! 


Data  Set  #1  includes  single-channel  waveform  data,  beams,  and  parametric  data 
from  241  regional  events  recorded  by  the  NORESS  and  ARCESS  arrays  in  Norway. 
The  purpose  for  assembling  this  data  set  is  to  use  it  to  develop  and  train  neural  net¬ 
works  to  perform  seismic  data  processing  and  interpretation  tasks.  However,  this  data 
set  may  also  be  useful  for  many  other  seismic  research  applications.  Data  Set  #1  is 
available  at  SAIC,  and  the  purpose  of  this  report  is  to  describe  it  in  detail. 

2.1  NORESS/ARCESS  Arrays 

The  NORESS  and  ARCESS  arrays  in  Norway  include  25  short-period  instruments 
in  four  concentric  rings  with  a  maximum  diameter  of  3  km  (Figure  1).  The  array 
configuration  and  sampling  rate  were  designed  to  enhance  the  detection  of  regional  sig¬ 
nals  [Mykkeltveit,  et  al.,  1983;  Mykkeltveit  and  Ringdal ,  1988].  The  radius  of  the  inner 
ring  (called  the  A-ring)  is  about  150  m.  The  radii  of  the  B-,  C-,  and  D-rings  are  300 
m,  700  m,  and  1500  m,  respectively.  The  number  of  sensors  on  the  A-,  B-,  C-,  and 
D-rings  are  3,  5,  7,  and  9,  respectively.  The  individual  station  locations  for  the 
NORESS  and  ARCESS  arrays  are  given  in  Table  2  (locations  are  given  relative  to  the 
reference  locations  listed  at  the  bottom  of  this  table).  Four  of  the  25  array  elements 
are  equipped  with  3-component  seismometers.  These  are  the  center  element  (A0),  and 
three  sensors  on  the  C-ring  (C2,  C4,  and  C7).  The  rest  of  the  array  elements  only 
have  vertical-component  seismometers. 

The  NORESS  and  ARCESS  data  are  continuously  recorded,  and  the  short-period 
data  are  digitized  at  a  rate  of  40  samples/s.  Figure  2  shows  the  short-period  instru¬ 
ment  response.  This  response  applies  to  all  elements  of  the  NORESS  and  ARCESS 
arrays.  The  instrument  response  is  approximately  fiat  to  velocity  between  2  and  8  Hz. 
The  digitization  gain  is  105  digital  counts/volt. 

2.2  Regional  Events 

This  section  describes  the  events  in  Data  Set  #1.  It  gives  the  event  locations, 
magnitudes,  epicentral  distances  from  NORESS  and  ARCESS,  and  identifications  (e.g., 
explosion  or  earthquake). 

2.2.1  Selection  Criteria 

Data  Set  #1  was  developed  to  support  the  seismic  application  of  DARPA’S 
Neural  Network  Program.  The  goal  of  this  application  is  determine  whether  or  not 
neural  networks  can  improve  upon  current  methods  for  seismic  monitoring  of  nuclear 
explosion  testing.  The  emphasis  is  on  low  yields,  so  the  primary  interest  in  on 
regional  distances  (<  20°).  We  used  data  from  NORESS  and  ARCESS  for  Data  Set 
#1  because  these  are  prototype  arrays  for  regional  monitoring,  and  there  is  a  large 
database  of  waveforms  and  parametric  data  archived  at  CSS.  The  NORESS  and 
ARCESS  data  are  continuously  processed  by  IMS,  and  the  results  of  the  automated 
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Figure  1.  The  location  and  array  geometry  are  plotted  for  the  NORESS  and 
ARCESS  arrays  (Figure  provided  by  Frode  Ringdal,  NORSAR). 


Table  2.  Station  Locations 


Station 

NORESS 

Elevation  dnorth1 
(km)  (km) 

deast1 

(km) 

Station 

ARCESS 

Elevation  dnorth2 
(km)  (km) 

deast2 

(km) 

NRAO 

.3020 

.0030 

.0040 

ARA0 

.4030 

.0010 

-.0003 

NRA1 

.2910 

.1460 

.0490 

ARA1 

.4110 

.1600 

.0530 

NRA2 

.3110 

-.1030 

.1080 

ARA2 

.3920 

-.1210 

.0770 

NRA3 

.2960 

-.0300 

-.1430 

ARA3 

.4020 

-.0300 

-.1490 

NRB1 

.2990 

.3210 

.0700 

ARBI 

.4140 

.3360 

.0820 

NRB2 

.3150 

.0300 

.3340 

ARB2 

.3970 

.0970 

.2940 

NRB3 

.3140 

-.2980 

.1430 

ARB3 

.3760 

-.2690 

.1890 

NRB4 

.2990 

-.2170 

-.2280 

ARB  4 

.3780 

-.2250 

-.2310 

NRB5 

.2890 

.1630 

-.2720 

ARB  5 

.4050 

.1580 

-.2830 

NRC1 

.2990 

.6870 

.1090 

ARC1 

.3810 

.6900 

.0810 

NRC2 

.3390 

.3410 

.6030 

ARC2 

.3950 

.3863 

.6657 

NRC3 

.3520 

-.2380 

.6470 

ARC3 

.3760 

-.2140 

.6730 

NRC4 

.3110 

-.6570 

.2080 

ARC4 

.3770 

-.6167 

.2287 

NRC5 

.2990 

-.5690 

-.3960 

ARC5 

.3740 

-.5380 

-.2960 

NRC6 

.3030 

-.0480 

-.6870 

ARC6 

.3950 

-.0810 

-.6830 

NRC7 

.2750 

.5480 

-.4470 

ARC7 

.3620 

.5300 

-.4700 

NRD1 

.3050 

1.4800 

.1920 

ARD1 

.3950 

1.4910 

.1350 

NRD2 

.3720 

1.0150 

1.0980 

ARD2 

.3660 

1.1430 

.9720 

NRD3 

.4530 

.0760 

1.4930 

ARD3 

.3310 

.1880 

1.6510 

NRD4 

.3790 

-.9010 

1.1890 

ARD4 

.3710 

-.8580 

1.1810 

NRD5 

.3480 

-1.4510 

.3350 

ARD5 

.3510 

-1.4940 

.2330 

NRD6 

.3520 

-1.3260 

-.6810 

ARD6 

.4130 

-1.3470 

-.6130 

NRD7 

.3370 

-.5660 

-1.3680 

ARD7 

.4130 

-.6070 

-1.3600 

NRD8 

.3010 

.4140 

-1.3360 

ARD8 

.3680 

.3920 

-1.4430 

NRD9 

.2780 

1.2570 

-.8020 

ARD9 

.3590 

1.1730 

-.7780 

1.  Relative  to  the  reference  location:  60.735°N,  11.541°E. 

2.  Relative  to  the  reference  location:  69.535°N,  25.506°E. 
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system  are  reviewed  by  seismic  analysts.  The  automated  results,  the  analyst  results, 
and  comparisons  between  them  are  stored  in  an  on-line  relational  database  at  CSS. 
Thus,  these  data  provide  an  excellent  opportunity  to  identify  deficiencies  in  the  current 
automated  system,  and  to  improve  its  performance. 

MIT  Lincoln  Laboratory  defined  several  goals  for  improving  the  automated  moni¬ 
toring  system  using  neural  network  techniques.  These  include  regional  phase 
identification,  onset  time  estimation,  typical  and  atypical  event  detection,  and  event 
identification  [LaCoss  et  al.,  1990;  LaCoss  et  at.,  1991].  Based  on  these  goals,  they 
requested  eight  separate  waveform  databases  for  Data  Set  #1.  These  are  listed  as 
DB1-DB3  in  Table  3.  They  also  requested  a  separate  parametric  database  (without 
waveforms)  for  10  days  of  continuous  operation.  This  is  listed  as  DB9  in  Table  3. 


Table  3.  Data  Set  #1 


Database 

Description 

Number  of  Events 
Requested  Provided 

DB1 

High-Quality  Regional  Analyst-Corrected 

50 

50 

DB2 

High-Quality  Regional  Analyst-Accepted 

50 

20 

DB3 

Random  Selection  of  Analyst-Corrected 

50 

50 

DB4 

Random  Selection  of  Analyst-Accepted 

50 

N/A 

DB5 

Random  Selection  of  Analyst-Rejected 

25 

N/A 

DB6 

Non-event  (Noise)  Detections 

40 

N/A 

DB7 

Teleseisms 

25 

N/A 

DB8 

Unusual  Events 

10 

21 

DB9 

10-days  of  Parametric  Data 

N/A 

249t 

DB10 

High-Quality  Regional  Analyst-Corrected 

N/A 

50 

DB1 1 

High-Quality  Regional  Analyst-Corrected 

N/A 

50 

t  Waveforms  were  not  provided. 


Data  Set  #1  was  developed  on  the  basis  of  this  request.  It  includes  waveforms 
and  parametric  data  from  the  Intelligent  Array  System  (IAS)  before  and  after  analyst 
review  (these  data  are  described  in  the  next  section).  IAS  was  the  predecessor  to  IMS. 
It  was  specifically  designed  for  the  two-array  network  of  NORESS  and  ARCESS,  and 
it  did  not  attempt  to  locate  teleseismic  events.  IMS  currently  processes  data  from  a 
network  of  four  European  arrays,  and  two  3-component  stations  in  Poland  (3- 
component  stations  in  the  Soviet  Union  and  China  will  also  be  added  to  this  network). 
IMS  locates  both  regional  and  teleseismic  events,  and  its  results  are  similar  to  those 
from  IAS  for  events  that  are  at  regional  distance  from  either  NORESS  or  ARCESS. 

We  used  data  from  IAS  rather  than  the  current  IMS  because  IMS  was  not  opera¬ 
tional  until  November  1990,  which  was  two  months  after  the  start  of  our  contract. 
Also,  the  rules  used  in  its  knowledge-based  system  were  modified  frequently  to 
improve  performance  until  approximately  April,  1991.  However,  results  from  IAS 
between  October,  1989  and  October,  1990  are  archived  at  CSS.  The  only  disadvantage 
of  using  the  IAS  data  is  that  waveforms  were  not  saved  for  all  events.  The  criteria  for 
saving  waveforms  varied  slightly  over  time,  but  they  were  based  primarily  on  the 
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event  location,  number  of  detecting  stations,  and  event  magnitude.  Because 
waveforms  were  saved  for  only  selected  events,  we  were  not  able  to  satisfy  all  the 
requests  in  DB1-DB8.  Therefore,  we  added  two  additional  databases,  DB10  and 
DB1 1. 

The  criteria  for  selecting  events  for  Data  Set  #1  varied  for  each  database  in  Table 
3.  For  all  databases,  we  required  event  latitudes  to  be  between  50°  and  80°,  and  longi¬ 
tudes  to  be  between  -10°  and  40°  (e.g.,  regional  distance  from  NORESS  or  ARCESS). 
We  also  required  that  the  waveforms  be  saved  on-line  at  CSS.  For  DB1,  DB10,  and 
DB11  we  require  "high-quality"  events  for  which  the  results  of  the  automated  system 
were  corrected  in  some  way  by  an  analyst  (e.g.,  retiming  a  detection,  renaming  a 
phase,  etc).  We  define  "high-quality"  as  events  detected  at  both  arrays  whose  locations 
are  constrained  by  at  least  3  defining  phases.  In  addition,  for  DB1  we  required  that  all 
events  be  detected  after  January  1,  1990  to  ensure  that  they  were  reviewed  by  the 
same  (NORSAR)  analyst.  This  constraint  had  to  be  relaxed  to  get  enough  events  for 
DB10  and  DB11.  The  same  criteria  were  used  for  DB3  (random  selection  of  analyst- 
corrected  events)  that  were  used  for  DB1,  except  there  were  no  constraints  on  the 
number  of  detecting  stations  or  the  number  of  defining  phases  (e.g.,  the  "high-quality" 
constraints). 

The  original  request  for  DB2  was  50  high-quality  events  for  which  the  results  of 
the  automated  system  were  accepted  by  an  analyst.  There  are  very  few  of  these  events 
with  waveforms  saved  in  the  IAS  database,  so  we  were  forced  to  relax  the  high-quality 
constraint.  Instead,  we  provided  all  analyst-accepted  events  with  waveforms  within 
our  latitude  and  longitude  bounds  (20  events).  Of  course,  this  meant  that  there  were 
no  events  to  provide  as  DB4.  The  request  for  DB5  was  for  events  that  were  formed 
by  the  automated  system,  and  then  rejected  by  an  analyst.  There  are  rejected  events  in 
the  IAS  database,  but  waveforms  are  not  saved.  Therefore,  we  could  not  provide  DB5. 
DB6  and  DB7  are  requests  for  noise  and  teleseismic  detections,  respectively.  A  detec¬ 
tion  is  labeled  "N"  for  noise  by  IAS  if  the  phase  velocity  estimated  using  a  broadband 
/  -k  method  is  less  than  2.8  km/s,  and  it  is  not  associated  witn  an  event.  Similarly,  a 
detection  is  labeled  "T"  for  teleseism  by  IAS  if  the  estimated  phase  velocity  is  greater 
than  14  km/s,  and  it  is  not  associated  with  a  regional  event.  Using  these  definitions, 
there  are  50  "teleseismic"  detections  and  142  "noise"  detections  in  the  other  waveform 
databases  in  Table  3.  However,  these  labels  are  based  only  on  the  estimate  of  phase 
velocity,  and  they  have  not  been  validated  by  an  analyst. 

The  request  for  DB8  was  for  10  unusual  events.  For  this,  we  gathered  data  from 
multiple-event  sequences  that  include  two  or  more  events  with  origin  times  that  are 
within  60  s  of  each  other.  These  events  produce  mixed  signals  (e.g..  interleaved 
phases)  at  either  NORESS  or  ARCESS,  or  both.  This  database  consists  of  21  events 
in  10  multiple-event  sequences.  The  only  other  criteria  applied  to  the  selection  of 
these  data  were  the  latitude  and  longitude  bounds. 

2.2.2  Event  Description 

Table  4  lists  the  241  events  in  Data  Set  #1  for  which  waveforms  were  provided. 
The  first  column  lists  the  orid,  which  is  a  unique  positive  integer  that  identifies  each 
event  in  the  CSS  database  [Anderson  et  al.,  1990].  The  second  column  lists  the 
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Table  4.  Events  With  Waveforms  in  Data  Set  #1 


orid 

Database 

Origin  Time 

191531 

DB1 

90  01  24 

16:21:47.014 

191546 

DB1 

90  01  25 

05:09:37.277 

191548 

DB1 

90  01  25 

09:49:22.456 

191552 

DB1 

90  01  25 

12:02:30.661 

191571 

DB1 

90  01  26 

10:01:19.878 

191576 

DB1 

90  01  26 

12:04:55.867 

191764 

DB1 

90  02  02 

08:42:06.695 

192008 

DB1 

90  02  08 

14:31:39.935 

191945 

DB1 

90  02  09 

09:32:50.929 

191954 

DB1 

90  02  09 

12:27:59.740 

192093 

DB1 

90  02  14 

10:16:05.013 

192387 

DB1 

90  02  23 

11:16:04.823 

192390 

DB1 

90  02  23 

12:05:06.429 

192671 

DB1 

90  03  05 

12:42:08.275 

192763 

DB1 

90  03  12 

14:34:31.648 

193041 

DB1 

90  03  23 

12:01:52.778 

193078 

DB1 

90  03  24 

15:28:54.658 

193185 

DB1 

90  03  28 

22:21:05.640 

193279 

DB1 

90  04  02 

13:41:10.653 

197468 

DB1 

90  04  09 

08:52:31.104 

197474 

DB1 

90  04  09 

13:02:30.367 

197517 

DB1 

90  04  11 

10:21:17.110 

197522 

DB1 

90  04  11 

11:00:23.168 

197535 

DB1 

90  04  11 

12:48:42.043 

197538 

DB1 

90  04  11 

13:46:06.790 

197576 

DB1 

90  04  12 

11:03:37.894 

197579 

DB1 

90  04  12 

12:09:06.463 

197628 

DB1 

90  04  14 

10:35:18.945 

197712 

DB1 

90  04  19 

11:26:55.068 

197742 

DB1 

90  04  20 

08:57:45.452 

197762 

DB1 

90  04  20 

12:23:40.040 

197764 

DB1 

90  04  20 

13:29:55.864 

197765 

DB1 

90  04  20 

13:56:21.070 

197868 

DB1 

90  05  02 

10:29:15.125 

198004 

DB1 

90  05  09 

16:25:44.411 

198020 

DB1 

90  05  10 

11:31:02.907 

198023 

DB1 

90  05  10 

11:53:46.825 

198052 

DB1 

90  05  11 

17:06:24.382 

198053 

DB1 

90  05  11 

17:07:24.053 

198143 

DB1 

90  05  15 

13:14:50.488 

198261 

DB1 

90  05  19 

20:40:12.834 

198271 

DB1 

90  05  20 

10:23:37.742 

198303 

DB1 

90  05  21 

11:30:47.309 

198384 

DB1 

90  05  22 

09:39:22.879 

198447 

DB1 

90  05  23 

13:37:47.510 

198461 

DB1 

90  05  23 

18:55:26.452 

198462 

DB1 

90  05  23 

19:09:30.078 

198523 

DB1 

90  05  26 

10:06:59.051 

198526 

DB1 

90  05  26 

11:13:00.712 

198546 

DB1 

90  05  27 

21:49:30.059 

103776 

DB2 

89  10  03 

12:29:34.655 

Latitude  Longitude 


3 

quake 

3 

quake  (H) 

4 

blast  (H) 

5 

blast  (H) 

6 

blast  (H) 

5 

blast  (H) 

7 

blast  (H) 

6 

blast  (H) 

4 

blast  (H) 

4 

blast 

6 

blast  (H) 

6 

blast  (H) 

5 

blast  (H) 

5 

expl 

5 

blast  (H) 

6 

blast  (H) 
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orid 

Database 

Origin  Time 

105012 

DB2 

89  10  04 

12:20:54.732 

105364 

DB2 

89  10  04 

19:19:44.234 

117198 

DB2 

89  10  12 

11:36:51.780 

118906 

DB2 

89  10  13 

11:25:12.799 

177958 

DB2 

89  11  21 

12:54:42.373 

184121 

DB2 

89  11  24 

13:01:27.445 

191018 

DB2 

90  01  07 

06:25:03.393 

191049 

DB2 

90  01  11 

10:10:01.214 

191654 

DB2 

90  02  02 

12:58:12.990 

191749 

DB2 

90  02  07 

07:54:16.686 

192034 

DB2 

90  02  15 

09:48:04.238 

192660 

DB2 

90  03  11 

14:14:32.868 

192764 

DB2 

90  03  14 

13:06:13.586 

193038 

DB2 

90  03  23 

10:12:38.438 

193094 

DB2 

90  03  27 

13:27:03.441 

197366 

DB2 

90  04  07 

05:23:43.882 

198311 

DB2 

90  05  23 

19:04:24.003 

200143 

DB2 

90  06  01 

12:01:24.068 

200149 

DB2 

90  06  01 

18:30:10.651 

191520 

DB3 

90  01  24 

12:01:54.714 

191526 

DB3 

90  01  24 

12:56:48.450 

191527 

DB3 

90  01  24 

13:31:56.883 

191583 

DB3 

90  01  26 

16:22:03.849 

191670 

DB3 

90  01  30 

12:02:20.415 

191806 

DB3 

90  02  04 

01:38:32.574 

191911 

DB3 

90  02  08 

11:31:13.041 

191955 

DB3 

90  02  09 

12:34:47.835 

191987 

DB3 

90  02  13 

11:24:56.765 

192192 

DB3 

90  02  16 

11:34:27.670 

192193 

DB3 

90  02  16 

11:38:47.495 

192199 

DB3 

90  02  16 

12:38:15.632 

192202 

DB3 

90  02  16 

12:55:04.285 

192380 

DB3 

90  02  23 

09:31:22.220 

192436 

DB3 

90  02  24 

19:37:18.179 

192482 

DB3 

90  02  26 

20:30:15.658 

192535 

DB3 

90  02  28 

13:30:42.498 

192554 

DB3 

90  03  01 

09:47:28.040 

192585 

DB3 

90  03  02 

11:01:56.012 

192587 

DB3 

90  03  02 

12:00:23.017 

192452 

DB3 

90  03  06 

10:47:54.178 

192454 

DB3 

90  03  06 

11:32:58.664 

192692 

DB3 

90  03  07 

10:24:51.491 

192791 

DB3 

90  03  14 

08:51:32.961 

192864 

DB3 

90  03  15 

11:54:23.370 

192887 

DB3 

90  03  16 

10:33:27.409 

192926 

DB3 

90  03  20 

13:36:59.973 

192969 

DB3 

90  03  21 

14:15:54.343 

193024 

DB3 

90  03  22 

13:21:51.699 

193034 

DB3 

90  03  23 

08:30:23.057 

193156 

DB3 

90  03  28 

10:59:13.988 

193166 

DB3 

90  03  28 

12:00:05.599 

193 187 

DB3 

90  03  29 

04:11:18.667 

193201 

DB3 

90  03  29 

11:27:18.825 

Latitude  Longitude  ML  nsta  ndef  Event  Type 


2  blast  (H) 

2  blast 

3  blast  (H) 

2  quake 

2  expl 

2  expl 

'  last 
last 
last 


last 

uak 


expl 

blast  (H) 
quake 
quake  (H) 
blast  (H) 
blast  (H) 
blast  (H) 
blast 
blast  (H) 
blast  (H) 
quake 
blast  (H) 
quake 
blast  (H) 
blast  (H) 
blast  (H) 
blast  (H) 
blast  (H) 
blast  (H) 
blast  (H) 
blast  (H) 
expl 

quake  (H) 
blast 
blast  (H) 
blast  (H) 
blast 
blast  (H) 
blast  (H) 
blast  (H) 

blast  (H) 
blast  (H) 
blast  (H) 
blast  (H) 
blast 
blast  (H) 
blast  (H) 
blast  (H) 
quake 
blast  (H) 
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orid 

Database 

Origin  Time 

193203 

DB3 

90  03  29 

11:31:43.507 

193218 

DB3 

90  03  30 

08:49:01.642 

193222 

DB3 

90  03  30 

10:39:03.617 

193250 

DB3 

90  03  31 

14:15:03.805 

193255 

DB3 

90  04  01 

05:06:31.114 

193266 

DB3 

90  04  02 

09:57:29.106 

193270 

DB3 

90  04  02 

10:31:55.798 

193280 

DB3 

90  04  02 

13:46:23.256 

193281 

DB3 

90  04  02 

14:03:34.589 

197381 

DB3 

90  04  03 

14:16:33.314 

197382 

DB3 

90  04  03 

14:24:33.159 

197408 

DB3 

90  04  06 

09:45:33.118 

197412 

DB3 

90  04  06 

10:36:01.684 

197413 

DB3 

90  04  06 

10:37:02.280 

197417 

DB3 

90  04  06 

11:34:29.834 

197432 

DB3 

90  04  06 

12:20:36.430 

106759 

DB8 

89  10  04 

10:12:40.564 

106762 

DB8 

89  10  04 

10:13:40.029 

118427 

DB8 

89  10  11 

06:22:29.539 

118499 

DB8 

89  10  11 

06:23:18.152 

131912 

DB8 

89  10  20 

14:55:52.442 

131916 

DB8 

89  10  20 

14:55:54.385 

131917 

DB8 

89  10  20 

14:56:01.685 

152005 

DB8 

89  11  02 

23:33:18.495 

152009 

DB8 

89  11  02 

23:33:23.372 

189277 

DB8 

89  11  21 

11:03:29.729 

189330 

DB8 

89  1 1  21 

11:03:57.490 

197533 

DB8 

90  04  11 

12:38:07.117 

197534 

DB8 

90  04  11 

12:39:15.714 

197593 

DB8 

90  04  13 

10:18:56.028 

197594 

DB8 

90  04  13 

10:19:41.850 

197601 

DB8 

90  04  13 

12:58:01.935 

197602 

DB8 

90  04  13 

12:59:18.645 

197704 

DB8 

90  04  19 

09:58:53.558 

197706 

DB8 

90  04  19 

10:00:42.828 

197818 

DB8 

90  04  28 

11:01:07.889 

197819 

DB8 

90  04  28 

11:02:11.118 

103539 

DB10 

89  10  02 

15:11:30.828 

118708 

DB10 

89  10  11 

07:49:36.455 

119051 

DB10 

89  10  11 

10:03:05.858 

123686 

DB10 

89  10  12 

12:20:51.001 

123713 

DB10 

89  10  12 

14:29:28.606 

123721 

DB10 

89  10  13 

11:41:30.794 

123728 

DB10 

89  10  13 

12:00:25.894 

125589 

DB10 

89  10  16 

10:21:15.181 

125900 

DB10 

89  10  16 

12:44:02.741 

127401 

DB10 

89  10  17 

15:00:47.687 

128464 

DB10 

89  10  17 

20:42:08.973 

131862 

DB10 

89  10  20 

06:41:23.940 

131911 

DB10 

89  10  20 

13:25:32.576 

134678 

DB10 

89  10  24 

15:10:35.164 

134688 

DB10 

89  10  24 

18:05:03.191 

135510 

DB10 

89  10  26 

09:24:12.089 

Latitude  Longitude 
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orid 

Database 

Origin  Time 

Latitude 

Longitude 

Ml 

nsta1 

nder 

Event  Type 

135514 

DB10 

89  10  26 

10:13:57.858 

61.19 

29.97 

246 

2 

6 

blast  (H) 

140720 

DB10 

89  10  26 

13:50:58.974 

60.54 

29.63 

2.42 

2 

6 

blast  (H) 

140725 

DB10 

89  10  26 

14:38:21.700 

61.56 

31.54 

2.42 

2 

6 

blast  (H) 

140759 

DB10 

89  10  27 

11:32:57.358 

64.81 

30.43 

2.36 

2 

6 

blast  (H) 

140775 

DB10 

89  10  27 

12:50:43.185 

69.45 

30.88 

2.65 

2 

6 

blast  (H) 

140819 

DB10 

89  10  28 

11:42:52.921 

61.85 

36.33 

2.58 

2 

6 

blast 

147829 

DB10 

89  11  01 

12:25:41.306 

59.16 

28.09 

2.66 

2 

6 

blast  (H) 

147845 

DB10 

89  11  01 

13:02:59.129 

60.93 

28.98 

2.57 

2 

6 

blast  (H) 

147846 

DB10 

89  1 1  01 

13:26:10.027 

60.88 

29.34 

2.68 

2 

5 

blast 

147858 

DB10 

89  11  01 

15:59:34  987 

61.90 

30.79 

2.54 

2 

6 

blast  (H) 

152019 

DB10 

89  11  03 

08:19:34  973 

67.68 

34.32 

2.87 

2 

7 

blast  (H) 

152059 

DB10 

89  11  03 

13:20:49.061 

61.60 

25.29 

2.35 

2 

5 

blast 

152097 

DB10 

89  11  04 

07:39:08.699 

67.89 

32.30 

2.28 

2 

5 

blast 

152109 

DB10 

89  11  04 

10:05:21.280 

64.78 

30.44 

2.70 

2 

6 

blast  (H) 

152944 

DB10 

89  11  04 

1 2:46:30.90  i 

69.36 

30.97 

3.03 

2 

5 

blast  (H) 

1529% 

DB10 

89  11  05 

11:05:54.398 

59.45 

27.11 

2.22 

2 

6 

blast  (H) 

153019 

DB10 

89  11  06 

05:58:40.963 

67.68 

33.56 

2.77 

2 

7 

blast  (H) 

153033 

DB10 

89  11  06 

11:59:14.428 

63.09 

28.24 

2.26 

2 

6 

blast  (H) 

155511 

DB10 

89  11  08 

22:56:30.837 

67.19 

20.% 

2.27 

2 

5 

blast  (H) 

160154 

DB10 

89  11  09 

13:28:18.765 

59.32 

25.21 

2.30 

2 

5 

blast  (H) 

161012 

DB10 

89  11  09 

14:12:21.951 

60.59 

29.17 

2.49 

2 

5 

blast  (H) 

160940 

DB10 

89  11  10 

07:06:14.404 

65.03 

26.41 

2.02 

2 

6 

quake  (H) 

160944 

DB10 

89  11  10 

10:42:14.172 

59.43 

27.27 

2.38 

2 

6 

blast  (H) 

161037 

DB10 

89  1 1  10 

12:03:11.747 

69.43 

30.93 

2.25 

2 

4 

blast  (H) 

160949 

DB10 

89  11  10 

12:40:38.302 

54.96 

16.40 

2.46 

2 

5 

expl 

161041 

DB10 

89  11  10 

13:43:38.951 

59.58 

25.36 

2.12 

2 

4 

blast  (H) 

161317 

DB10 

89  11  11 

07:10:27.217 

67.63 

33.56 

2.17 

2 

6 

blast  (H) 

161049 

DB10 

89  11  11 

08:41:27.120 

67.64 

33.97 

2.52 

2 

7 

blast  (H) 

172081 

DB10 

89  11  14 

12:00:33.595 

61.17 

30.06 

2.29 

2 

4 

blast  (H) 

172095 

DB10 

89  11  14 

14:47:22.230 

59.48 

26.50 

2.76 

2 

6 

blast  (H) 

172363 

DB10 

89  1 1  15 

12:08:12.985 

63.29 

27.74 

2.15 

2 

6 

blast  (H) 

175063 

DB10 

89  1 1  16 

12:35:06.559 

58.93 

18.71 

2.04 

2 

5 

expl 

175305 

DB10 

89  1 1  16 

14:55:49.832 

59.53 

27.14 

2.24 

2 

5 

blast  (H) 

177371 

DB10 

89  11  17 

11:00:01.966 

64.80 

30.72 

2.08 

2 

6 

blast  (H) 

107824 

DB1 1 

89  09  29 

08:23:16.393 

67.53 

33.58 

2.61 

2 

5 

blast  (H) 

107849 

DB11 

89  09  29 

12:02:26.738 

60.97 

29.31 

2.87 

2 

6 

blast  (H) 

107883 

DB1 1 

89  09  29 

12:30:30.804 

61.32 

28.77 

2.42 

2 

6 

blast  (H) 

107888 

DB1 1 

89  09  29 

12:54:34.980 

69.40 

30.86 

2.26 

2 

3 

blast  (H) 

107890 

DB11 

89  09  29 

13:00:00.756 

59.47 

26.80 

2.54 

2 

6 

blast  (H) 

107914 

DB11 

89  09  29 

14:13:34.956 

59.53 

24.% 

2.30 

2 

4 

blast  (H) 

108561 

DB11 

89  09  30 

13:10:08.696 

67.63 

29.84 

2.49 

2 

4 

blast  (H) 

105219 

DB11 

89  10  03 

11:56:42.915 

64.16 

23.67 

2.15 

2 

4 

expl 

106770 

DB1 1 

89  10  04 

1 1:27:41 .636 

61.58 

30.33 

2.29 

2 

5 

blast  (H) 

106780 

DB11 

89  10  04 

13:16:20.806 

55.79 

16.75 

2.30 

2 

4 

exp) 

107868 

DB1 1 

89  10  05 

11:49:15.124 

59.46 

23.86 

2.17 

2 

4 

blast 

108213 

DB1 1 

89  10  05 

14:48:27.340 

59.59 

26.47 

2.54 

2 

6 

blast  (H) 

111577 

DB1 1 

89  10  06 

09:05:55.448 

67.59 

35.01 

2.71 

2 

4 

blast  (H) 

111585 

DB1 1 

89  10  06 

10:16:39.153 

64.77 

31.20 

2.44 

2 

5 

blast  (H) 

111787 

DB1 1 

89  10  06 

1 1: 12:39.704 

59.43 

27.13 

2.71 

2 

6 

blast  (H) 

111946 

DB1 1 

89  10  06 

11:42:14.650 

67.64 

29.85 

2.26 

2 

3 

blast  (H) 

112156 

DB11 

89  10  06 

12:01:19.638 

63.24 

27.59 

2.12 

2 

7 

blast  (H) 

115620 

DB11 

89  10  09 

16:20:14.188 

67.08 

21.24 

2.02 

2 

4 

blast  (H) 

116470 

DB1 1 

89  10  10 

13:25:58.927 

59.02 

25.52 

2.81 

2 

5 

blast  (H) 
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orid 

Database 

Origin  Time 

Latitude 

Longitude 

Ml 

nsta1 

ndef2 

Event  Type 

190865 

DB11 

89  10  17 

14:55:04.120 

59.55 

26.78 

2.56 

2 

5 

blast  (H) 

133590 

DB11 

89  10  23 

13:28:47.537 

61.85 

31.12 

2.44 

2 

6 

blast  (H) 

140713 

DB11 

89  10  26 

12:59:05.936 

59.67 

28.22 

2.26 

2 

5 

blast  (H) 

140826 

DB11 

89  10  27 

15:10:30.505 

67.72 

33.83 

2.12 

2 

4 

blast  (H) 

190881 

DB1I 

89  11  16 

13:17:56.963 

2.33 

2 

4 

blast  (H) 

177383 

DB11 

89  11  17 

08:24:54.337 

67.67 

34.15 

2.75 

2 

7 

blast  (H) 

177438 

DB11 

89  11  18 

11:26:46.268 

61.86 

2.59 

2 

6 

blast  (H) 

178839 

DB11 

89  11  20 

12:57:57.410 

2.24 

2 

5 

blast  (H) 

189367 

DB11 

89  11  21 

13:05:42.979 

61.27 

29.78 

2.28 

2 

6 

blast  (H) 

189479 

DB11 

89  11  24 

10:00:34.632 

64.84 

2.33 

2 

4 

blast  (H) 

189461 

DB11 

89  11  24 

12:04:08.162 

67.63 

2.91 

2 

7 

blast 

189469 

DB11 

89  11  24 

16:02:09.780 

66.95 

21.66 

2.35 

2 

5 

blast  (H) 

189655 

DB11 

89  11  25 

11:11:45.657 

33.52 

2.27 

2 

7 

blast  (H) 

189658 

DB11 

89  11  25 

12:37:40.038 

67.62 

2.37 

2 

5 

blast  (H) 

192440 

DB11 

90  02  25 

10:00:26.090 

6.83 

2.29 

2 

5 

quake  (H) 

197524 

DB11 

90  04  11 

11:13:08.723 

59.17 

27.37 

2.29 

2 

6 

blast  (H) 

197580 

DB1 1 

90  04  12 

12:12:15.582 

28.98 

2.18 

2 

5 

blast  (H) 

197596 

DB1 1 

90  04  13 

10:47:07.274 

28.33 

2.33 

2 

4 

blast 

197639 

DB11 

90  04  14 

15:15:22.056 

34.43 

WEM 

2 

expl 

197671 

DB11 

90  04  18 

10:12:54.031 

59.29 

27.89 

Ha 

2 

6 

blast  (H) 

197674 

DB11 

90  04  18 

11:39:32.142 

27.74 

2.39 

2 

4 

blast  (H) 

197700 

DB11 

90  04  19 

08:05:27.018 

69.16 

34.64 

2.45 

2 

6 

blast  (H) 

197741 

DB11 

90  04  20 

08:38:51.774 

67.61 

33.73 

2.51 

2 

7 

blast  (H) 

197747 

DB11 

90  04  20 

09:57:37.358 

27.16 

2.32 

2 

6 

blast  (H) 

197820 

DB1 1 

90  04  28 

12:02:02.118 

67.98 

33.76 

2 

5 

blast  (H) 

197825 

DB11 

90  04  28 

14:37:06.044 

59.51 

26.46 

2.23 

2 

6 

blast  (H) 

197876 

DB1 1 

90  05  02 

12:05:59.365 

57.74 

11.72 

3.83 

2 

4 

expl 

198017 

DB11 

90  05  10 

11:02:46.221 

59.33 

27.22 

2.29 

2 

6 

blast  (H) 

198039 

DB1 1 

90  05  11 

11:11:46.391 

59.29 

2.33 

2 

5 

blast 

198272 

DB1 1 

90  05  20 

10:27:07.366 

2 

5 

quake  (H) 

198353 

DB1 1 

90  05  21 

12:50:39.845 

28.11 

2.39 

2 

4 

blast 

1.  The  number  of  detecting  stations. 

2.  The  number  of  defining  phases  (number  of  phases  used  to  locate  the  event). 
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database  name  from  Table  3.  The  rest  of  the  columns  list  the  event  origin  time,  lati¬ 
tude,  longitude,  Lg  magnitude  (ML),  number  of  detecting  stations,  number  of  defining 
phases,  and  event  type.  Events  that  were  identified  in  a  regional  bulletin  produced  by 
the  University  of  Helsinki  are  indicated  by  "(H)"  following  the  event  type  in  Table  4. 
Other  events  were  identified  by  Sereno  [1991]. 

The  event  locations  are  plotted  in  Figure  3.  Of  the  241  events,  215  were 
recorded  at  both  NORESS  and  ARCESS,  and  221  had  at  least  3  defining  phases.  The 
epicentral  distances  are  primarily  between  200  and  1800  km  (Figure  4).  The  Lg  mag¬ 
nitude  distribution  is  plotted  in  Figure  5.  This  magnitude  is  computed  from  the  peak 
amplitude  on  a  2-4  Hz  incoherent  beam  in  the  time  window  defined  by  group  veloci¬ 
ties  of  3.0  and  3.6  km/s  [Bache  et  ai,  1991].  The  magnitudes  of  the  events  in  Data 
Set  #1  are  primarily  between  1.5  and  3.2. 

Most  of  the  events  in  Data  Set  #1  were  identified  (mine  blast  or  earthquake)  in 
the  regional  bulletin  produced  by  the  University  of  Helsinki.  The  other  events  were 
either  not  identified  by  the  University  of  Helsinki,  cm-  they  were  not  in  their  bulletin. 
In  either  case,  these  events  were  identified  by  Sereno  [1991]  on  the  basis  of  location, 
origin  time,  S/P  amplitude  ratios,  spectral  variance,  and  past  seismicity.  Along  with 
the  identification  of  each  of  these  events,  Sereno  [1991]  gives  a  brief  description  of  the 
basis  for  the  identification  and  a  subjective  measure  of  confidence  (high,  medium,  or 
low).  Table  5  lists  the  number  of  earthquakes,  mine  blasts,  underwater  explosions,  and 
other  explosions  for  each  of  the  waveform  databases  in  Table  3.  Events  that  are  prob¬ 
able  explosions  and  are  located  onshore  (but  not  near  known  mines)  are  labeled  "other 
explosions." 


Table  5.  Event  Identification  for  Data  Set  #1 


Number  of  Events 

Event  Type 

DB1 

DB2 

DB3 

DB8 

DB10 

DB11 

Total 

Earthquakes 

5 

4 

5 

0 

2 

2 

18 

Mine  Blasts 

34 

11 

39 

10 

43 

44 

181 

Underwater  Explosions 

8 

1 

3 

4 

3 

3 

22 

Other  Explosions 

3 

3 

2 

7 

2 

1 

18 

Not  Identified 

0 

1 

1 

0 

0 

0 

2 

Total 

50 

20 

50 

21 

50 

50 

241 

Figure  6  shows  the  locations  of  the  presumed  earthquakes,  mine  blasts,  and  other 
explosions  (either  underwater  or  onshore).  Most  of  the  earthquakes  are  located  on  the 
Mid-Atlantic  Ridge,  or  near  the  west  coast  of  Norway.  The  mine  blasts  are  primarily 
in  the  northwestern  USSR  (e.g.,  Estonia,  Leningrad,  Kola  Peninsula),  northern  Sweden, 
and  Finland.  Most  of  the  underwater  explosions  are  located  in  the  Baltic  Sea. 

The  locations  of  the  249  events  in  DB9  are  shown  in  Figure  7.  These  events 
were  recorded  over  a  continuous  10-day  period  starting  February  28,  1990.  Data  Set 
#1  does  not  include  waveform  data  for  these  events,  and  the  parametric  data  arc 
described  in  the  next  section. 
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Events  With  Waveforms  in  Data  Set  #1 


ARCESS 


Distance  (Km) 

4.  Histograms  of  epicentra)  distance  to  the  events  in  Data  Set  #1  are  plotted 
ARCESS,  and  (b)  NORESS. 


Magnitude  Distribution  -  Data  Set  #1 

St - 


0.0  0.5  1.0  13  2.0  23  3.0  33  4.0 

IMS  Magnitude  (ML) 


Figure  5.  A  histogram  of  Lg  magnitude  is  shown  few  events  in  Data  Set  #1.  The 
mean  ML  is  2.31,  and  the  standard  deviation  is  0.42. 
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Event  Identification  -  Data  Set  #1 


°  Earthquake  ■  Mine  Blast  •  Explosion 


Figure  6.  This  map  shows  the  locations  of  earthquakes  (circles),  mine  blasts 
(squares),  and  other  explosions  (asterisks)  in  Data  Set  #1. 
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Events  With  Only  Parametric  Data  (DB9)  -  Data  Set  #1 


Figure  7.  Epicenters  are  plotted  for  the  249  events  with  only  parametric  data  (DB9) 
in  Data  Set  #1. 
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2.3  Data  Exchange  Format 


This  section  describes  the  exchange  format  for  Data  Set  #1.  These  data  were  pro¬ 
vided  to  MU  Lincoln  Laboratory  on  read/write  optical  disks,  but  they  are  also  avail¬ 
able  on  9-track  tapes  in  UNIX  tar  format.  Database  reports  that  include  hard-copy 
displays  and  information  for  each  event  were  also  provided  to  MIT  Lincoln  Labora¬ 
tory.  For  each  event,  these  reports  include  a  map  showing  station  and  event  location, 
a  listing  of  event  origin  and  detection  data,  the  standard  waveform  display  that  is  pro¬ 
vided  to  the  IAS  seismic  analyst,  and  plots  of  each  detecting  and  3-component  beam 
(described  below).  Separate  reports  were  provided  for  DB1,  DB2,  DB3,  DB8,  DB9, 
DB10,  and  DB11  [Patnaik  and  Sereno,  1990a]. 

Data  Set  #1  includes  parametric  data  generated  by  the  automated  IAS  processing, 
and  all  changes  made  by  seismic  analysts.  It  also  includes  all  available  short-period 
waveform  data  recorded  at  NORESS  and  ARCESS  when  the  same  event  is  detected  by 
both  arrays.  The  parametric  data  are  stored  as  ASCII  files,  and  the  waveform  data  are 
stored  in  the  binary  format  used  by  SAC  ( Seismic  Analysis  Code).  SAC  is  a  widely- 
used  interactive  analysis  software  package  that  was  developed  by  Joseph  Tull  at 
Lawrence  Livermore  National  Laboratory.  Section  2.3.1  briefly  describes  the 
parametric  data.  The  formats  of  the  ASCII  data  files  are  described  in  detail  in  Appen¬ 
dix  A.  The  waveform  data  are  described  in  Section  2.3.2.  An  example  of  the  displays 
and  data  that  are  included  in  each  database  report  is  given  in  Section  2.3.3. 

2.3.1  Parametric  Data 

Figure  8  shows  the  directory  structure  for  Data  Set  #1  (directories  are  printed  in 
bold-face  type,  and  files  are  printed  in  italics).  The  top-level  directory  includes  sub¬ 
directories  for  each  of  the  databases  listed  in  Table  3,  and  sub  directories  for  documen¬ 
tation  and  static  data  (e.g.,  data  that  are  the  same  for  all  events).  The  documentation 
directory  includes  a  description  of  the  data  exchange  format,  and  a  description  of 
several  minor  enhancements  that  we  made  to  SAC.  The  static  data  include  beam 
recipe  files  (described  below),  the  NORESS/ARCESS  short-period  instrument  response, 
the  travel-time  tables  for  regional  phases  used  in  IAS,  and  the  locations  of  the  indivi¬ 
dual  stations  in  each  array. 

The  general  directory  structure  for  each  database  is  shown  in  the  middle  of  Figure 
8.  Each  database  directory  includes  up  to  four  sub-directories,  and  two  ASCII  files. 
The  README  file  gives  information  for  a  few  specific  events.  The 
ExpSys  Analyst.dbX  file  summarizes  corrections  made  by  the  analyst  to  the  results  of 
the  automated  processing  (the  variable  X  appended  to  a  file  name  or  directory  stands 
for  the  database  number).  For  example,  it  includes  the  distance  between  the  locations 
determined  by  the  expert  system  and  the  analyst,  the  difference  in  their  origin  times, 
and  the  number  of  phases  that  were  added  or  retimed  by  the  analyst  (see  Appendix  A). 

The  ExpSys  directory  contains  the  results  of  the  automated  IAS  processing.  This 
includes  the  four  files  in  Figure  8  with  the  prefix  IEB,  which  stands  for  initial  event 
bulletin.  lEB.orig  lists  event  origin  data  (latitude,  longitude,  depth,  and  origin  time). 
lEB.det  lists  detection  and  association  data  such  as  phase  identification,  arrival  time. 


20 


DIRECTORY  STRUCTURE  /  DATA  SET  #1 


Top-Level  Directory 

DB1  DB2  DB3  DB8  DB9  DB10  DB11  doc  static 


Database  Directory 


DBX 

README  ExpSys 

ExpSys  Analyst. dbX  IEB.orig 

lEB.det 

lEB.apma 

lEBsbsnr 


Analyst 

FEB.orig 

FEB.det 

FEB.distaz 


SAC  Waveform  Directory 


SAC/oridXXXXX 

ARAO.se.sac 
ARAO.sn.sac 
ARAO.sz.sac 
ARAl.sz.sac 
ARA2szsac 
ARA3.sz.sac 
ARBI  sz  sac 
ARB2.sz.sac 
ARB3.sz.sac 
ARB4.sz.sac 
ARB5.sz.sac 
ARCl.sz.sac 
ARC2sesac 
ARC 2. sn. sac 
ARC2.sz.sac 
ARC3.sz.sac 


ARC4.se.sac 
ARC4.sn.sac 
ARC4.sz.sac 
ARC5.sz.sac 
ARC6.szsac 
ARC7.se.sac 
ARC7.sn.sac 
ARC7.sz.sac 
ARDlsz.sac 
ARD2.sz.sac 
ARD3sz.sac 
ARD4sz.sac 
ARDS  sz. sac 
ARD6.sz.sac 
ARD7sz.sac 
ARD8.sz.sac 


ARD9.szsac 

NRAOsesac 

NRAO.snsac 

NRAOsz.sac 

NRAl.sz.sac 

NRA2.sz.sac 

NRA3sz.sac 

NRBlsz.sac 

NRB2sz.sac 

NRB3sz.sac 

NRB4sz.sac 

NRBSsz.sac 

NRC1  .szsac 

NRC2.se  sac 

NRC2.sn.sac 

NRC2.szsac 


Figure  8.  The  directory  structure  is  described 


EVID 
README 
EVID.dbX 
Helsinki  orig 
Helsinkibul.dbX 
MSMP.dbX 
CEPPKS.dbX 
SPVAR.dbX 


SAC 

oridXXXXX 

oridXXXXX 

oridXXXXX 

etc 


NRC3sz.sac 
NRC4sesac 
NRC4snsac 
NRC4sz.sac 
NRC5sz.sac 
NRC6sz.sac 
NRC7sesac 
NRC7  snsac 
NRC7sz.sac 
NRDlszsac 
NRD2szsac 
NRD3szsac 
NRD4szsac 
NRD5  rzsac 
NRD6szsac 
NRD7szsac 


for  Data  Set  #1. 


NRD8sz.sac 

NRD9sz.sac 

ARC.cbxxx.sac 

ARC.cbxxx.sac 

ARC.cbxxxsac 

ARC.cbxxx  sac 

ARC.cbxxxsac 

ARC.cbxxxsac 

ARC.cbxxxsac 

NOR.cbxxxsac 

NOR.cbxxxsac 

NOR.cbxxxsac 

NOR.cbxxxsac 

NOR.cbxxxsac 
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azimuth  and  phase  velocity  estimated  from  f  -k  processing,  amplitude,  signal-to-noise 
ratio,  and  frequency.  lEB.apma  lists  results  from  automated  particle  motion  analysis 
[Jurkevics,  1988].  This  includes  estimates  of  rectilinearity,  planarity,  long-  and  short- 
axis  incidence  angles,  and  horizontal-to-vertical  power  ratio  for  each  detection. 
lEB.sbsnr  gives  signal  and  noise  amplitudes  measured  on  a  standard  set  of  six  beams. 

The  Analyst  directory  contains  parametric  data  after  review  by  a  seismic  analyst. 
The  files  FEB.orig  and  FFB.det  are  similar  to  the  corresponding  files  under  ExpSys, 
but  they  include  changes  made  by  the  analyst  (the  FEB  prefix  stands  for  final  event 
bulletin).  FEB.distaz  lists  the  epicentral  distance  and  station-to-everit  azimuth  from 
NORESS  and  ARCESS  for  each  event  in  FEB.orig.  The  Analyst  directory  does  not 
include  particle  motion  or  standard-beam  amplitude  files  since  these  attributes  are  not 
recalculated  after  analyst  review. 

The  EVID  directory  contains  the  identification  (e.g.,  earthquake  or  explosion)  of 
each  event  in  FEB.orig.  This  identification  is  bat^d  primarily  on  a  regional  seismic 
bulletin  produced  by  the  University  of  Helsinki.  The  README  file  gives  information 
for  a  few  specific  events.  EVID.dbX  gives  the  identification  of  each  event, 
Helsinki.orig  gives  the  origin  information  from  the  Helsinki  Bulletin,  and 
Helsinkibul.dbX  gives  the  complete  unedited  listing  from  Helsinki  Bulletin.  The  other 
three  files  in  this  directory  contain  data  that  are  relevant  for  event  identification. 
MSMP.dbX  lists  regional  P- wave  magnitudes  computed  from  the  amplitude  of  Pn  an 
Pg ,  and  regional  S  -wave  magnitudes  computed  from  the  amplitude  of  Sn  and  Lg .  The 
difference  between  these  magnitudes  is  a  possible  discriminant  (high  values  of  ms  -  mp 
indicate  the  event  is  an  earthquake,  low  values  are  inconclusive).  The  two  files  called 
CEPPKS.dbX  and  SPVAR.dbX  give  the  results  of  cepstral  analysis,  and  are  useful  for 
identifying  ripple-fired  mining  explosions  | Baumgardt  and  Ziegler ,  1987],  These  data 
are  described  in  Appendix  A. 

2.3.2  Waveform  Data 

The  SAC  sub  directory  contains  the  NORESS  and  ARCESS  waveform  data  in 
SAC  format.  There  is  a  separate  directory  for  each  event.  These  directories  are 
labeled  as  oridXXXXXX,  where  XXX XXX  refers  to  the  unique  integer  origin 
identification  in  FEB  orig.  There  are  separate  SAC  data  files  for  each  short-period 
channel  in  the  NORESS  and  ARCESS  arrays  (33  channels/array).  These  waveform 
data  files  are  7-minute  segments  that  start  30  seconds  before  the  theoretical  Pn  arrival 
time  (based  on  the  final  event  origin).  The  files  are  named  as  station. channel. sac 
where  station  is  the  station  code  for  each  array  element  (e.g.,  NRAO,  NRA1,  NRA2, 
etc),  and  channel  is  the  channel  code  (sz  is  short-period  vertical,  se  is  short-period 
east,  and  sn  is  short-period  north). 

In  addition  to  the  single-channel  waveforms,  each  directory  includes  several 
coherent  beams  in  SAC  format.  Coherent  beams  are  formed  by  steering  the  single- 
channel  waveforms  using  a  specified  velocity  and  azimuth,  stacking,  and  filtering  over 
a  specified  frequency  band.  These  files  are  named  as  array. cbxxx. sac,  where  array  is 
either  NCR  for  NORESS  or  ARC  for  ARCESS,  cb  stands  for  coherent  beam,  and  xxx 
is  the  beam  number  in  the  first  column  of  Table  6  (this  table  describes  the  beam  set 
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Table  6.  IAS  Beam  Deployment 
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Beam 

Velocity 
(k m/s) 

Filter 

(Hz) 

Filter 

Order 

Azimuth 

(degrees) 

Beam 

Type1 

Ring  Subset 

289 

OO 

5.0-10.0 

3 

0.0 

C 

A0 

B 

C 

5.0-10.0 

3 

30.0 

C 

A0 

B 

C 

291 

5.0-10.0 

3 

90.0 

C 

A0 

B 

C 

292 

5.0-10.0 

3 

150.0 

C 

A0 

B 

C 

293 

5.0-10.0 

3 

C 

A0 

B 

C 

294 

5.0-10.0 

3 

C 

A0 

B 

C 

295 

5.0-10.0 

3 

330.0 

C 

A0 

B 

C 

296 

9.9 

8.0-16.0 

3 

30.0 

c 

A0 

A 

B 

297 

9.9 

8.0-16.0 

3 

90.0 

c 

A0 

A 

B 

298 

9.9 

8.0-16.0 

3 

150.0 

c 

A0 

A 

B 

299 

9.9 

8.0-16.0 

3 

210.0 

c 

A0 

A 

B 

9.9 

8.0-16.0 

3 

270.0 

c 

AO 

A 

B 

301 

9.9 

8.0-16.0 

3 

330.0 

c 

AO 

A 

B 

15.9 

1. 5-3.5 

3 

80.0 

c 

AO 

C 

D 

303 

15.9 

3 

80.0 

c 

AO 

C 

D 

304 

15.9 

2.5-4.5 

3 

80.0 

c 

AO 

B 

C 

D 

305 

15.9 

3 

80.0 

c 

AO 

B 

C 

D 

1. 5-3.5 

3 

30.0 

c 

AO 

C 

D 

2.0-4 .0 

3 

30.0 

c 

AO 

C 

D 

308 

2.5-4.5 

3 

30.0 

c 

AO 

B 

C 

D 

309 

3 

30.0 

c 

AO 

B 

C 

D 

310 

OO 

2 

0.0 

I 

AO 

C 

312 

OO 

3 

0.0 

I 

AO 

C 

313 

OO 

2 

0.0 

I 

AO 

C 

1.  C  =  Coherent,  I  =  Incoherent  (vertical  channels),  H  =  Incoherent  (horizontal  channels) 
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applied  to  NORESS  and  ARCESS  during  the  operation  of  IAS).  The  detecting  beam  is 
defined  as  the  standard  beam  in  Table  6  with  the  highest  snr.  This  beam  is  included 
for  each  detection  that  is  associated  with  a  final  origin.  However,  if  the  detecting 
beam  is  incoherent,  then  it  is  replaced  by  a  coherent  beam  that  uses  the  same  beam- 
forming  parameters. 

Coherent  beams  that  are  calculated  using  data  from  the  four  3-component  ele¬ 
ments  of  each  array  are  also  included  in  Data  Set  #1.  Three  of  these  beams  are  calcu¬ 
lated  for  each  detection  that  is  associated  with  a  final  event  origin  (one  for  each  com¬ 
ponent;  vertical,  north-south,  and  east-west).  These  3-component  beams  use  the  same 
beamforming  parameters  as  the  detecting  beam,  except  that  the  array  subset  only 
includes  the  four  3-component  elements.  The  number  assigned  to  the  beam  formed 
from  the  vertical  components  is  200  plus  the  beam  number  of  the  detecting  beam. 
Similarly,  the  numbers  assigned  to  the  beams  formed  from  the  north- south  and  east- 
west  components  are  400  and  600  plus  the  beam  number  of  the  detecting  beam, 
respectively. 

2.33  Database  Reports 

This  section  gives  an  example  of  the  displays  and  information  that  are  included  in 
the  database  reports  provided  to  MIT  Lincoln  Laboratory  [Patnaik  and  Sereno,  1991a]. 
This  example  is  for  a  mining  explosion  in  western  USSR.  Figure  9  is  a  map  showing 
the  locations  of  NORESS  and  ARCESS,  and  the  event  epicenter.  Table  7  lists  detec¬ 
tion  data  before  analyst  review  (e.g.,  the  results  of  the  automated  IAS  processing),  and 
Table  8  lists  detection  data  after  analyst  review.  Figure  10  shows  the  dislay  beams  for 
ARCESS  (this  is  the  standard  waveform  display  that  is  provided  to  the  MS  seismic 
analyst).  Figure  1 1  plots  the  detecting  beam  for  Pn  at  ARCESS,  and  the  3-component 
beams  (vertical,  north-south,  and  east- west).  Similarly,  Figure  12  plots  the  beams  for 
Sn  at  ARCESS.  Beams  are  not  included  for  Lg  since  this  phase  was  added  by  the 
analyst  (e.g.,  it  was  not  detected  by  the  automated  processing).  Figure  13  shows  the 
display  beams  for  NORESS,  and  Figure  14  shows  the  detecting  and  3-component 
beams  for  Pn  at  NORESS.  Beams  were  not  computed  for  Sn  and  Lg  phases  at 
NORESS  because  they  were  added  by  the  analyst. 
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ORID=  192093  90/02/14  10  16  5.013  LAT  61.70  LON  31.37  ML  2.2 


Figure  9.  This  map  shows  the  NORESS  and  ARCESS  station  locations,  and  the 
epicenter  of  an  event  in  Date  Set  #1  (after  analyst  review). 
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ARID  ST*  CHAM  CRANIO  XI  MM  DO  HR:  Ml:  S3  .US  IPHASC  PHASK 
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(5)  Add  an  Sn  phase  at  NORESS. 

(6)  Add  an  Lg  phase  at  NORESS. 

The  location  determined  by  the  expert  system  is  about  250  km  from  the  by  location 
determined  by  the  analyst  for  this  event. 


ORID=  192093  ARCESS 
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ORID  192093,  Coherent  Beams:  detecting,  vert, 
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Figure  11.  This  figure  shows  the  detecting  beam  for  Pn  at  ARCESS  in  the  top 
panel,  and  the  beams  formed  from  the  3 -component  array  elements  in  the  bottom 
panels  (vertical,  north-south,  east-west). 


ORID  192093,  Coherent  Beams:  detecting,  vert, 
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Figure  12.  This  figure  shows  the  detecting  beam  for  Sn  at  ARCESS  in  the  top 
panel,  and  the  beams  formed  from  the  3-component  array  elements  in  the  bottom 
panels  (vertical,  north-south,  east-west). 


ORID=  192093  NORESS 


NORESS  display  beams  arc  plotted  for  the  event  in  Figure  9.  Associ 
nd  Lg  phases  are  highlighted  (see  caption  for  Figure  10). 


ORID  192093,  Coherent  Beoms:  detecting,  vert 
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Figure  14.  This  figure  shows  the  detecting  beam  for  Pn  at  NORESS  in  the  top 
panel,  and  the  beams  formed  from  the  3-component  array  elements  in  the  bottom 
panels  (vertical,  north-south,  east-west). 


(THIS  PAGE  INTENTIONALLY  LEFT  BLANK) 
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3.  SUMMARY 

Wc  developed  a  data  set  to  test  and  evaluate  the  performance  of  neural 
networks  for  automated  processing  and  interpretation  of  seismic  data.  This  data  set 
may  also  be  valuable  for  many  other  studies  related  to  seismic  monitoring  of  nuclear 
explosion  testing.  It  includes  waveform  and  parametric  data  from  241  regional 
events  recorded  by  the  short-period  elements  of  the  NORESS  and  ARCESS  arrays  in 
Norway  (33  channels/array).  The  waveform  data  are  stored  in  SAC  binary  format, 
and  the  parametric  data  are  stored  in  ASCII  files.  The  event  epicentral  distances  are 
200-1800  km,  and  the  event  Lg  magnitudes  are  approximately  1. 5-3.2.  Most  of  the 
events  are  mining  explosions  in  western  USSR,  Sweden,  and  Finland.  However,  18  of 
the  events  are  earthquakes,  and  22  are  presumed  underwater  explosions.  Detailed 
documentation  has  been  developed  for  each  event,  and  is  included  in  eight  separate 
database  reports.  The  data  and  software  are  available  at  the  Center  for  Seismic 
Studies. 
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APPENDIX  A:  PARAMETRIC  DATA  EXCHANGE  FORMAT 


This  appendix  describes  the  parametric  data  files  in  Data  Set  #1  (Figure  2.8). 
These  data  were  retrieved  from  the  IAS  relational  database  at  the  Center  for  Seismic 
Studies.  Examples  of  each  data  file  are  given  at  the  end  of  this  section. 

ExpSys  Analyst.dbX 


Table  A.l  ExpSys _Analyst.dbX  Data  File 


ExpSys_AnaIyst.dbX 

attribute 

field 

Storage 

attribute 

name 

no. 

type 

description 

eorid 

1 

i4 

expert  system  origin  id 

forid 

2 

i4 

final  origin  id 

ddist 

3 

f4 

distance  between  forid  and  eorid 

dtime 

4 

f4 

origin  time  difference 

rprimp 

5 

i4 

no.  of  renamed  primary  phases 

rsecondp 

6 

i4 

no.  of  renamed  secondary  phases 

added 

7 

i4 

no.  of  added  phases 

retime 

8 

i4 

no.  of  retimed  phases 

Table  A.l  describes  the  ExpSys  Analyst.dbX  data  file.  This  file  summarizes  the 
corrections  made  by  an  analyst  to  the  results  of  the  automated  processing.  All  of  the 
attributes  in  this  table  were  extracted  from  the  ex  an  database  table  at  CSS  which  is 
an  IMS  extension  to  the  Center  Version  3  Database  [Swanger  et  al.,  1991].  This  file 
joins  the  ExpSys  data  files  and  the  Analyst  data  files  through  the  unique  integer  origin 
identifications,  eorid  and  forid.  The  attributes  have  the  following  definitions: 


eorid: 
forid : 
ddist: 
dtime: 


rprimp: 
rsecondp : 
added: 

retime: 


A  unique  positive  integer  which  identifies  an  origin  determined  by  the 
expert  system. 

A  unique  positive  integer  which  identifies  an  origin  determined  (or  vali¬ 
dated)  by  a  seismic  analyst. 

The  distance  between  the  event  location  determined  by  the  expert  system 
and  the  event  location  determined  by  the  analyst  (in  kilometers). 

The  absolute  value  of  the  difference  between  the  origin  time  estimated  by 
the  expert  system  and  the  origin  time  estimated  by  the  analyst  (in 
seconds). 

Number  of  primary  (P  -type)  phases  that  were  renamed  by  an  analyst. 
Number  of  secondary  (5 -type)  phases  that  were  renamed  by  an  analyst 

Number  of  phases  that  were  added  by  an  analyst.  These  are  phases  that 
were  not  detected  by  the  automated  system. 

Number  of  phases  that  were  retimed  by  an  analyst. 
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There  are  cases  in  which  fields  5-8  are  zero,  but  the  analyst  and  expert  system 
solutions  differ.  This  can  happen  if  the  analyst  associates  or  disassociates  a  phase  with 
the  current  event  (the  phase  is  not  renamed,  retimed,  or  added). 

ExpSys 

The  ExpSys  directory  contains  parametric  data  from  the  results  of  the  automated 
processing  (before  analyst  review).  The  files  in  this  directory  have  the  prefix  1EB  for 
initial  event  bulletin. 

IEB.orig 


Table  A.2  IEB.orig  Data  File 


IEB.orig 

attribute 

field 

Storage 

attribute 

name 

no. 

type 

description 

i4 

expert  system  origin  id 

date 

date  (yy  mm  dd) 

time 

time  (hnmm:ss.ms) 

f4 

estimated  latitude 

5 

f4 

estimated  longitude 

6 

f4 

estimated  depth 

Table  A.2  describes  the  IEB.orig  data  file.  This  file  gives  the  event  location  and 
origin  time  determined  by  the  expert  system.  The  link  between  the  expert  system  ori¬ 
gin  ( eorid)  and  die  final  analyst-reviewed  origin  (forid)  is  given  in  the 
ExpSys  Analyst. dbX  data  file.  All  of  the  attributes  in  this  file  were  extracted  from  die 
origin  Version  3  database  table  at  CSS  [Anderson  et  al„  1991],  The  attributes  have 
the  following  definitions: 

eorid :  A  unique  positive  integer  which  identifies  an  origin  determined  by  the 

expert  system. 

date:  Date  in  year-month-day  format.  This  is  converted  from  the  Julian  date 

format  in  the  CSS  database. 

time:  Event  origin  time  (GMT).  This  is  converted  from  the  epoch  time  in  the 

CSS  database. 

lot:  Estimated  event  latitude.  Latitudes  are  positive  north  of  the  equator. 

Ion :  Estimated  event  longitude.  Longitudes  are  positive  east  of  the  Greenwich 

meridian. 

depth:  Estimated  source  depth  in  kilometers.  The  depths  of  all  events  in  Data  Set 

#1  were  constrained  to  0.0  km. 

Magnitudes  are  not  calculated  for  the  solutions  determined  by  the  expert  system. 
They  are  calculated  after  analyst  review  (see  FEB.orig). 
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IEBdet 


Table  A. 3  IEB.det  Data  File 


IEB.det 

attribute 

field 

Storage 

attribute 

name 

no. 

■KE3I 

description 

eorid 

1 

i4 

expert  system  origin  id 

arid 

2 

i4 

arrival  id 

sta 

3 

c6 

station  code 

chan 

4 

c8 

channel  code 

chanid 

5 

i4 

channel  id 

date 

6 

date 

date  (yy  mm  dd) 

time 

7 

time 

time  (hnmm:ss.ms) 

iphase 

8 

c8 

reported  phase 

phase 

9 

c8 

associated  phase 

amp 

10 

f4 

amplitude  (nm) 

ffeq 

11 

f4 

frequency 

snr 

12 

f4 

signal-to-noise  ratio 

velo 

13 

f4 

observed  phase  velocity 

azimuth 

14 

f4 

observed  azimuth 

fkq 

15 

i4 

f-k  quality  measure 

Table  A.  3  describes  the  lEB.det  data  file.  This  file  gives  the  detection  and  phase 

association  data  determined  by  the  expert  system.  The  attributes  in  this  file  were 

extracted  from  the  arrival  and  assoc  Version  3  database  tables  at  CSS  [Anderson  et 

al.,  1991].  The  attributes  have  the  following  definitions: 

eorid:  A  unique  positive  integer  which  identifies  an  origin  determined  by  the 

expert  system. 

arid:  A  unique  positive  integer  which  identifies  an  arrival. 

sta:  Code  name  given  to  a  seismic  station.  Station  codes  for  NORESS  and 

ARCESS  are  NRAO  and  ARAO,  respectively. 

chan:  Code  name  given  to  the  data  channel.  Code  names  for  Data  Set  #1  are  zb 

and  hb  for  beams  formed  from  data  recorded  by  vertical-  and  horizontal- 
component  sensors,  respectively. 

chanid:  A  unique  positive  integer  which  identifies  the  recording  channel.  For  Data 

Set  #1  this  is  equal  to  the  IAS  beam  number  in  Table  2.5. 

date:  Date  in  year-month-day  format.  This  is  converted  from  the  Julian  date 

format  in  the  CSS  database. 

time:  Arrival  time  (GMT).  This  is  converted  from  the  epoch  time  in  the  CSS 

database. 
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iphase :  The  name  initially  given  to  a  seismic  phase.  In  Data  Set  #1  (and  in  the 

IAS  database  at  CSS)  this  is  equal  to  the  final  phase  identification  (phase , 
see  below)  if  it  is  associated  with  an  event.  Otherwise,  it  is  die  phase 
name  assigned  by  die  expert  system  during  single-station  processing.  This 
is  discussed  in  more  detail  below. 

phase:  Name  given  to  each  associated  phase  (e.g.,  final  phase  identification  after 

network  processing). 

amp:  This  is  zero- to- peak  displacement  amplitude  measured  on  the  detecting 

channel  in  nanometers.  It  is  corrected  for  the  instrument  response  at  the 
dominant  signal  frequency  (see  freq  below). 

freq:  Dominant  signal  frequency  in  Hertz.  This  is  equal  to  the  reciprocal  of  the 

period,  per ,  in  the  CSS  Version  3  Database  [Anderson  et  al.,  1990]. 

snr:  Signal-to-noise  ratio  measured  on  the  detecting  channel.  This  is  equal  to 

the  short-term-average  amplitude  (STAV)  divided  by  the  long-term- 
average  amplitude  (LTAV)  immediately  preceding  the  detection  [Bache  et 
al.,  1990]. 

veto:  Phase  velocity  estimated  using  a  broadband  f-k  method  (km/s).  This  is 

equal  to  the  reciprocal  of  the  slowness,  slow,  in  the  IAS  database  at  CSS. 

azimuth:  Azimuth  estimated  using  a  broadband  f-k  method  (measured  in  degrees 

clockwise  from  north). 

fkq:  Integer  measure  of  the  quality  of  the  f-k  spectrum.  Possible  values  are 

1-4,  with  the  highest  quality  indicated  by  fkq  -  1  [Bache  et  al.,  1990]. 

The  IEB.det  file  includes  all  detections  in  a  7-minute  window  starting  30  s  before 
the  theoretical  Pn  arrival  time  for  each  event  (based  on  the  final  origin  in  FEB.orig). 
This  includes  phases  that  are  associated  with  the  event,  those  that  are  associated  with  a 
different  event  (e.g„  mixed  events),  and  those  that  are  not  associated  with  any  event. 
If  the  detection  is  not  associated  with  an  event,  then  eorid  is  set  to  -1. 

The  possible  values  for  iphase  in  Data  Set  #1  are  Pn,  Pg,  Px,  P ,  Sn,  Lg,  Sx,  S, 
T,  and  iV.  Px  and  P  are  synonyms,  as  are  Sx  and  S.  These  labels  are  used  for  P  -  and 
S-type  phases  for  which  the  path  from  source  to  receiver  is  unknown  (e.g.,  coda  detec¬ 
tions).  The  label  T  is  assigned  to  detections  with  estimated  phase  velocities  >14  km/s 
(probable  teleseisms).  N  is  used  if  the  phase  velocity  is  <  2.8  km/s,  or  if  fkq  =  4 
(probable  "noise"  detections).  The  final  label  assigned  to  phases  that  are  associated 

with  an  event,  phase,  is  Pn,  Pg,  Px,  P ,  Sn,  Lg,  Sx,  or  S.  This  label  is  set  to  " - " 

for  detections  that  are  not  associated  with  a  regional  event. 

The  IAS  definitions  of  iphase  and  phase  are  slightly  different  than  the  correspond¬ 
ing  IMS  definitions  because  IMS  also  detects  and  locates  teleseismic  events.  For 
example,  P  stands  for  teleseismic  P  in  the  IMS  database,  and  it  is  not  synonymous 
with  Px.  Also,  iphase  is  not  equal  to  phase  for  associated  detections.  Instead,  it  is 
equal  to  the  phase  label  that  is  assigned  during  single-station  processing. 
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lEB.apma 

Table  A. 4  lEB.apma  Data  File 


IEItapma 

attribute 

field 

Storage  attribute 

name 

no. 

HE713 

description 

phase 

1 

c8 

phase 

arid 

2 

i4 

arrival  id 

freq 

3 

f4 

frequency 

snr 

4 

f4 

signal-to-noise  ratio 

amp  3 

5 

f4 

3-component  amplitude 

rect 

6 

f4 

rectilinearity 

plan 

7 

f4 

planarity 

hvrat 

8 

f4 

horizontal-to-vertical  ratio 

hmxmn 

9 

f4 

max-to-min  horizontal  ratio 

inang3 

10 

f4 

shon-axis  incidence  angle 

azimuth 

11 

f4 

observed  azimuth 

ema 

12 

f4 

emergence  angle 

prime 

13 

f8 

P  phase  extraction  time 

srime 

14 

f8 

S  phase  extraction  rim? 

auth 

15 

cl5 

author 

commid 

16 

i4 

comment  id 

Table  A.4  describes  the  IEB.apma  data  file.  This  file  gives  results  from 
automated  particle  motion  analysis.  All  attributes  in  this  file  were  extracted  from  the 
apma  database  table  at  CSS  which  is  an  IMS  extension  to  the  Center  Version  3  Data¬ 
base  [Swanger  et  al.,  1991]. 

The  method  used  in  IMS  for  particle  motion  analysis  was  developed  by  Jurkevics 
[1988].  It  computes  the  polarization  ellipse  within  overlapping  time  windows  by  solv¬ 
ing  the  eigenproblem  for  the  covariance  matrix.  Data  from  the  four  3-component  sen¬ 
sors  at  NORESS  and  ARCESS  are  combined  by  averaging  the  individual  covariance 
matrices  before  solving  the  eigenproblem.  The  covariance  matrices  are  computed  in 
the  time  domain  for  several  frequency  bands,  and  then  normalized  and  averaged  to 
obtain  a  wide-band  estimate  for  each  of  the  overlapping  windows.  The  IAS  implemen¬ 
tation  of  this  is  described  by  Bache  et  al.  [1990]. 

Several  of  the  particle  motion  attributes  are  calculated  from  the  time  window  with 
the  maximum  rcctilinearity.  These  are  called  P- type  attributes  in  the  following 
description.  Also,  several  attributes  are  calculated  from  the  time  window  with  max¬ 
imum  3-component  amplitude.  The  are  called  S-type  attributes.  The  panicle  motion 
attributes  have  the  following  definitions: 

phase:  Name  given  to  each  associated  phase.  This  is  currently  not  filled  in  the 

IAS  database,  so  it  is  set  to  the  N/A  value  of  for  Data  Set  #1. 

arid:  A  unique  positive  integer  which  identifies  an  arrival. 
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freq: 


snr. 


amp3: 


rect 


plan: 


hvrat. 


hmxmn: 


Center  frequency  of  the  passbands  with  snr  >  1.5  used  in  the  averaging. 
The  passbands  are  1-2,  2-4,  4-8  and  8-16  Hz.  For  example,  if  all  bands 
had  snr  >  1.5,  then  freq  would  be  8.5  Hz. 

Average  snr  for  frequency  bands  that  contribute  to  the  final  polarization 
estimates.  Each  snr  is  the  ratio  of  the  maximum  signal  3 -component 
amplitude  to  the  maximum  pre-arrival  noise  3-component  amplitude  (see 
amp3  below). 

Three-component  amplitude  measured  from  the  time  window  with  max¬ 
imum  rectilinearity  (P-type  attribute).  amp3  is  equal  to  the  sum  of  the 
square  roots  of  the  eigenvalues  (e.g.,  it  is  the  sum  of  the  amplitudes  meas¬ 
ured  along  the  3  axes  of  the  polarization  ellipsoid).  It  is  called  ampp  in 
later  versions  of  the  IMS  database. 

Signal  rectilinearity  defined  as: 


rect 


^•3  +  ^2 

2X, 


where  X,,  Xj,  and  Xj  are  the  eigenvalues  such  that  X,  >  X,  >  X,.  Later  ver¬ 
sions  of  IMS  use  the  square  root  of  this  quantity,  rect  is  a  A* -type  attri¬ 
bute. 

Signal  planarity  defined  as: 


plan 


Later  versions  of  IMS  use  the  square  root  of  this  quantity.  Planarity  is 
measured  from  the  window  with  the  maximum  3-component  amplitude 
(S-type  attribute).  It  is  called  plans  in  later  versions  of  the  IMS  database. 

Horizontal  to  vertical  power  ratio  defined  as: 


hvrat 


Cj  +  C2 

2c, 


where  c,,  c2,  and  c3  are  the  diagonal  elements  of  the  covariance  matrix, 
and  c,  corresponds  to  the  vertical  component,  hvrat  is  an  S-type  attribute. 

Maximum  to  minimum  horizontal  amplitude  ratio  defined  as: 

hmxmn  — 


where  X,  and  X2  are  the  maximum  and  minimum  eigenvalues  obtained  by 
solving  the  2-D  eigensystem  using  the  only  the  horizontal  components.  It 
is  an  S  -type  attribute. 
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inang3 : 

azimuth: 

ema: 

ptime : 
stime : 
auth: 
commid: 


Incidence  angle  (measured  from  the  vertical)  of  the  eigenvector  associated 
with  the  smallest  eigenvalue.  It  is  also  called  the  short-axis  incidence 
angle,  and  it  is  an  5 -type  attribute. 

Azimuth  of  the  eigenvector  associated  with  the  largest  eigenvalue.  It  is 
corrected  by  180°  to  give  an  estimate  of  the  station-to-event  azimuth. 
azimuth  is  a  P- type  attribute. 

Apparent  incidence  angle  (measured  from  the  vertical)  of  the  eigenvector 
associated  with  the  largest  eigenvalue.  It  is  also  called  the  long-axis 
incidence  angle,  or  the  emergence  angle.  It  is  called  inangl  in  later  ver¬ 
sions  of  the  IMS  database.  It  is  a  P  -type  attribute. 

Center  of  the  time  window  with  maximum  rectilinearity.  All  P- type  attri¬ 
butes  are  measured  from  this  time  window. 

Center  of  the  time  with  maximum  3-component  amplitude.  All  5 -type 
attributes  are  measured  from  this  time  window. 

Author.  This  field  is  set  to  the  N/A  value  of  for  all  arrivals  in  Data 
Set  #1. 

Comment  identification.  This  attribute  is  set  to  the  N/A  value  of  -1  for 
arrivals  in  Data  Set  #1. 


IEBsbsnr 


Table  A.5  IEBsbsnr  Data  File 


IEBjsbsnr 

attribute 

field 

Storage 

attribute 

name 

no. 

type 

description 

arid 

1 

i4 

arrival  id 

chanid 

2 

i4 

channel  id 

stav 

3 

f4 

short-term  average 

ltav 

4 

f4 

long-term  average 

Table  A.5  describes  the  IEBsbsnr  data  file.  This  file  gives  short-term  average 
signal  amplitudes  (STAV)  and  long-term  average  noise  amplitudes  (LTAV)  on  up  to 
six  standard  beams  (these  are  beams  201,  207,  254,  282,  310,  and  312  in  Table  2.5). 
All  attributes  in  this  file  were  extracted  from  die  sbsnr  database  table  at  CSS  which  is 
an  IMS  extension  to  the  Center  Version  3  Database  [Swanger  et  al.,  1991].  The  attri¬ 
butes  have  the  following  definitions: 

arid:  A  unique  positive  integer  which  identifies  an  arrival. 

chanid:  A  unique  positive  integer  which  identifies  the  recording  channel.  For  Data 

Set  #1  this  is  equal  to  the  IAS  beam  number  in  Table  2.5. 
stav:  Short-term  average  signal  amplitude  in  digital  counts  [Bache  et  al.,  1990]. 

ltav:  Long-term  average  pre-signal  amplitude  in  digital  counts  [Bache  et  al., 

1990]. 
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The  IEB.sbsnr  table  includes  the  amplitudes  on  all  six  beams  for  detections  that 
axe  associated  with  an  event.  However,  only  the  amplitudes  on  beam  312  are  saved 
for  unassociated  detections.  Also,  there  are  some  cases  where  all  beams  could  not  be 
formed  because  of  missing  or  bad  data. 

Analyst 


The  Analyst  directory  contains  parametric  data  after  review  by  a  seismic  analyst 
The  files  in  this  directory  have  the  prefix  FEB  for  final  event  bulletin. 


FEBx>rjg 


Table  A.6  FEB.orig  Data  Hie 


FEB.orig 

attribute 

field 

Storage 

attribute 

name 

no. 

type 

description 

forid 

1 

i4 

final  origin  id 

date 

2 

date 

date  (yy  mm  dd) 

time 

3 

time 

time  (hrrmmrss.ms) 

lat 

4 

f4 

estimated  latitude 

Ion 

5 

f4 

estimated  longitude 

depth 

6 

f4 

estimated  depth 

ml 

7 

f4 

local  magnitude 

nsta 

8 

i4 

No.  of  recording  stations 

ndef 

9 

i4 

No.  of  defining  phases 

Table  A.6  describes  the  FEB.orig  data  file.  This  file  gives  the  event  location  and 
origin  time  after  analyst  review  of  the  expert  system  solution.  The  attributes  in  this 
file  were  extracted  from  the  origin  Version  3  database  table  at  CSS  [Anderson  et  al., 
1991],  and  the  ev_summary  database  table  which  is  an  IMS  extension  to  the  Center 
Version  3  Database  [Swanger  et  al.,  1991]. 

The  first  6  attributes  in  this  file  have  the  same  definitions  as  those  described 
under  lEB.orig  except  that  the  expert  system  origin  identification  is  replaced  with  the 
final  origin  identification.  The  other  three  attributes  have  the  following  definitions: 

ml:  Local  Lg  magnitude  (see  below),  magnitude). 

nsta  Number  recording  stations. 

ndef  Number  of  defining  phases  (e.g.,  phases  that  are  used  to  constrain  the 
event  location). 

The  local  magnitude  is  referred  to  as  the  IMS  Version  1  MLg  by  Bache  et  al. 
[1991].  It  is  computed  from  the  peak  amplitude  on  a  2-4  Hz  incoherent  beam  in  the 
time  window  defined  by  group  velocities  of  3.0  to  3.6  km/s.  It  is  computed  for  each 
array  contributing  any  associated  phase  to  the  location  solution  (i.e.,  this  magnitude  is 
computed  even  when  there  is  no  detected  and  identified  Lg  phase). 


FEB.det 


Table  A.7  FEB.det  Data  File 


FEB.det 

attribute 

field 

Storage 

attribute 

name 

no. 

type 

description 

forid 

1 

i4 

final  origin  id 

arid 

2 

i4 

arrival  id 

sta 

3 

c6 

station  code 

chan 

4 

c8 

channel  code 

chanid 

5 

i4 

channel  id 

date 

6 

date 

date  (yy  mm  dd) 

time 

7 

time 

time  (hr.mmtss.ms) 

iphase 

8 

c8 

reported  phase 

phase 

9 

c8 

associated  phase 

amp 

10 

f4 

amplitude  (nm) 

freq 

11 

f4 

frequency 

snr 

12 

f4 

signal-to-noise  ratio 

velo 

13 

f4 

observed  phase  velocity 

azimuth 

14 

f4 

observed  azimuth 

fkq 

15 

i4 

f-k  quality  measure 

Table  A.7  describes  the  FEB.det  data  file.  This  file  gives  the  detection  and  phase 
association  data  after  analyst  review  of  the  automated  processing.  The  attributes  in 
this  file  were  extracted  from  the  arrival  and  assoc  Version  3  database  tables  at  CSS 
[Anderson  et  al .,  1991]. 

The  attributes  in  this  file  have  the  same  definitions  as  those  described  under 
lEB.det  except  that  the  expert  system  origin  identification  is  replaced  with  the  final  ori¬ 
gin  identification.  Phases  that  were  detected  by  the  automated  system  appear  in  both 
lEB.det  and  FEB.det  under  the  same  unique  arrival  identification,  arid.  The  only 
difference  between  the  two  files  for  these  detections  is  that  the  time  and/or  phase  label 
may  have  been  changed  by  the  analyst.  All  other  parameters  are  the  same  in  both 
files.  Phases  that  are  added  by  an  analyst  (that  were  not  detected  by  the  automated 
system)  appear  only  in  FEB.det.  These  phases  have  null  values  for  fields  4-5  and 
10-15  since  signal  processing  is  not  recalled. 
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FEBAistaz 


Tabic  A.8  FEB.distaz  Data  File 


FEB.distaz 

attribute 

field 

Storage 

attribute 

name 

no. 

HKTSSi 

description 

forid 

1 

i4 

final  origin  id 

sta 

2 

c6 

station  code 

distance 

3 

f4 

epicentral  distance  (km) 

seaz 

4 

f4 

station-to-event  azimuth 

Table  A.8  describes  the  FEB.distaz  data  file.  This  file  gives  the  distance  and 
station-to-event  azimuth  to  NORESS  and  ARCESS  for  each  event  in  FEB.orig.  These 
were  calculated  from  the  locations  of  the  center  elements  of  each  array,  and  the  event 
locations  in  FEB.orig. 

EVID 


The  EVID  directory  contains  the  identification  (earthquake  or  explosion)  of  each 
event  This  identification  is  based  primarily  on  a  regional  seismic  bulletin  produced  by 
the  University  of  Helsinki  [ Sereno ,  1991]. 


EVID.dbX 


Table  A.9  EVID.dbX  Data  File 


EVID.dbX 

attribute 

field 

Storage 

attribute 

name 

no. 

■KTsSHI 

description 

forid 

1 

i4 

final  origin  id 

evtype 

2 

c!5 

event  type 

Table  A.9  describes  the  EVID.dbX  data  file.  This  file  contains  the  probable 
identification  of  each  event  in  FEB.orig.  The  attributes  have  the  following  definitions: 

forid :  A  unique  positive  integer  which  identifies  an  origin  determined  (or  vali¬ 

dated)  by  a  seismic  analyst. 

evtype :  Event  type.  This  is  earthquake,  mine  blast,  or  explosion  (see  Section  2.2 

of  this  report).  If  the  event  type  is  followed  by  ”(H)",  then  the  event  was 
identified  in  the  Helsinki  Bulletin.  Otherwise,  the  event  was  identified  by 
Sereno  [1991]. 
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Helsinki.orig 


Table  A.  10  Helsinki.orig  Data  File 


Helsinki.orig 

attribute 

field 

Storage 

attribute 

name 

no. 

jgmnam 

description 

forid 

1 

i4 

final  origin  id 

date 

2 

date 

date  (yy  mm  dd) 

time 

3 

time 

time  (hr:mm:ss.ms> 

lat 

4 

f4 

estimated  latitude 

Ion 

5 

f4 

estimated  longitude 

depth 

6 

f4 

estimated  depth 

ml 

7 

f4 

local  magnitude 

evtype 

o 

O 

cl5 

event  type 

Table  A.  10  describes  the  Helsinki. orig  data  file.  This  file  contains  the  origin  data 
from  the  Helsinki  Bulletin  for  the  IMS  events  that  were  reported  in  that  bulletin.  The 
complete  unedited  listing  from  the  Helsinki  Bulletin  is  given  for  each  event  in 
Helsinkibul.dbX.  Most  of  the  attributes  in  this  table  were  described  previously.  The 
others  have  the  following  definitions: 

ml:  Local  magnitude  from  ~ie  Helsinki  Bulletin.  This  duration-based  magni¬ 

tude  is  described  by  Wahlstr'dm  and  Ahjos  [1984], 

evtype:  Event  type  determined  by  the  University  of  Helsinki.  This  is  either  man¬ 

ioc  for  manual  location  (these  are  mining  explosions),  earthquake,  or  an 
N/A  value  of  " - "  if  the  event  type  was  not  determined. 


MSMP.dbX 

Table  A.  11  MSMP.dbX  Data  File 


MSMP.dbX 

attribute 

field 

Storage 

attribute 

name 

no. 

type 

description 

forid 

1 

i4 

final  origin  id 

ml 

2 

f4 

local  magnitude 

mip 

3 

f4 

regional  P-wave  magnitude 

mis 

4 

f4 

regional  S-wave  magnitude 

msmp 

5 

f4 

mis  -  mip 

Table  A.  11  describes  the  MSMP.dbX  data  file.  This  file  lists  regional  P  -wave 
magnitudes  computed  from  Pn  and  Pg  amplitudes,  and  regional  S- wave  magnitudes 
computed  from  Sn  and  Lg  amplitudes.  The  difference  between  them,  msmp,  is  a  pos¬ 
sible  discriminant  (high  values  of  this  difference  indicate  that  the  event  is  an  earth- 


49 


quake,  and  low  values  are  inconclusive).  The  first  two  attributes  are  described  under 
FEB.orig.  The  others  have  the  following  definitions: 


mlp:  Regional  P- wave  magnitude  computed  from  Pn  and  Pg  amplitudes  [Bache 

et  al.y  1991;  IMS  Version  2  magnitudes]. 

mlp:  Regional  S  -wave  magnitude  computed  from  Sn  and  I.g  amplitudes  [Bache 

et  al.,  1991;  IMS  Version  2  magnitudes]. 
msmp :  The  difference  between  mis  and  mlp. 

CEPPKS.dbX 

Table  A.12  CEPPKSAbX  Data  File 


CEPPKS.dbX 


attribute 

name 

field 

no. 

Storage 

type 

attribute 

description 

forid- 

1 

i4 

final  origin  iu 

sta 

2 

c6 

station  code 

ptyp 

3 

c6 

cepstral  peak  type  code 

pkqf 

4 

f4 

cepstral  peak  quefrency 

Table  A.12  describes  the  CEPPKS.dbX  data  file.  This  file  gives  the  results  of 
cepstral  analysis  and  is  useful  for  identifying  ripple-fired  mining  explosions 
[Baumgardt  and  Zeigler,  1987].  The  first  two  attributes  were  described  previously. 
The  others  have  the  following  definitions: 

ptyp :  Consistent  cepstral  peak  type.  This  is  FC-PHS  if  consistent  Fourier  cep¬ 

stral  peaks  are  found  across  two  or  more  phases  for  one  array,  and  there  is 
no  peak  in  the  noise  cepstrum  at  this  quefrency.  Otherwise,  it  is  if  no 
consistent  cepstral  peaks  arc  found. 

pkqf:  Quefrency  of  the  consistent  cepstral  peak  (in  seconds).  This  is  set  to  zero 

if  there  are  no  consistent  peaks. 

The  best  evidence  for  ripple-firing  is  consistent  cepstral  peaks  across  phases  and 
arrays.  These  are  identified  in  the  CEPPKS.dbX  file  as  peaks  with  the  same  quefrency 
at  both  arrays.  Cepstral  peaks  that  appear  in  only  one  phase  at  a  given  station  are  not 
reported  as  consistent  peaks,  even  if  there  is  no  peak  in  the  noise  cepstrum  at  that 
quefrency. 


SPVAR.dbX 


Table  A.  13  SPVARxibX  Data  File 


SPVAR.dbX 

attribute 

field 

Storage 

attribute 

name 

no. 

description 

forid 

1 

i4 

final  origin  id 

arid 

2 

i4 

arrival  id 

sta 

3 

c6 

station  code 

phase 

4 

c8 

associated  phase 

acoef 

5 

f4 

"a"  coefficient  for  non-linear  trend 

bcoef 

6 

f4 

"b"  coefficient  for  non-linear  trend 

ccoef 

7 

f4 

"c"  coefficient  for  non-linear  trend 

fmin 

8 

f4 

min  frequency 

fmax 

9 

f4 

max  frequency 

svar 

10 

f4 

variance  of  detrended  log  spectrum 

Table  A.  13  describes  the  SPVAR.dbX  data  file.  This  file  gives  the  variance  of  the 
detrended  log  spectrum  for  each  phase  that  is  associated  with  an  event.  The  first  four 
attributes  were  described  previously.  The  others  have  the  following  definitions: 

( a,b,c)coef :  Three  coefficients  of  the  quadratic  trend  of  the  log  spectrum  between  fre¬ 
quencies  fmin  and  fmax. 

fmin\  Minimum  frequency  of  a  band  with  snr  >  3  dB. 

fmax:  Maximum  frequency  of  a  band  with  snr  >  3  dB. 

svar:  Variance  of  the  detrended  log  spectrum  between  fmin  and  fmax. 

Array-averaged  spectra  are  computed  for  a  5-s  window  staring  0.3  s  before  each 
arrival.  Noise  spectra  are  calculated  for  a  5-s  window  starting  5.3  s  before  the  first 
arrival.  All  spectra  are  corrected  for  the  short-period  NORESS/ARCESS  instrument 
response.  A  log  snr  spectrum  is  calculated  for  each  detection  using  the  noise  before 
the  first  arrival.  A  running  mean  (width  =  0.75  Hz)  is  applied  to  this  snr  spectrum, 
and  frequency  bands  fmin  to  fmax  are  determined  such  that  the  smoothed  snr  spectrum 
is  >  3  db  (only  frequency  bandwidths  >  4  Hz  are  retained).  A  second-order  polyno¬ 
mial  is  fit  to  the  log  signal  spectra  between  fmin  and  fmax  (the  coefficients  are  acoef, 
bcoef,  and  ccoef),  and  this  trend  is  removed.  The  variance  of  the  detrended  log  spec¬ 
trum  is  calculated  and  written  to  the  table  as  svar. 

Example  Data  Files 

An  example  of  each  parametric  data  file  is  given  in  the  next  few  pages  for  a  min¬ 
ing  explosion  in  western  USSR. 
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